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ABSTRACT 

Diffuse radio emission in the form of radio halos and relics has been found in a number of 
merging galaxy clusters. These structures indicate that shock and turbulence associated with 
the merger accelerate electrons to relativistic energies. We report the discovery of a radio relic 
+ radio halo system in PSZ1 G108.18-11.53 (z = 0.335). This cluster hosts the second most 
powerful double radio relic system ever discovered. We observed PSZ1 G108.18-11.53 with 
the Giant Meterwave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Tele¬ 
scope (WSRT). We obtained radio maps at 147, 323, 607 and 1380 MHz. We also observed 
the cluster with the Keck telescope, obtaining the spectroscopic redshift for 42 cluster mem¬ 
bers. From the injection index we obtained the Mach number of the shocks generating the two 
radio relics. For the southern shock we found M = 2.33^ ^, while the northern shock Mach 
number goes from M = 2.20^ ^ in the north part down to M = 2.00^ ^ in the southern re¬ 
gion. If the relation between the injection index and the Mach number predicted by diffusive 
shock acceleration (DSA) theory holds, this is the first observational evidence for a gradient 
in the Mach number along a galaxy cluster merger shock. 

Key words: galaxies: clusters: individual: PSZ1 G108.18-11.53 - large-scale structure of 
Universe - radio continuum: general 


1 INTRODUCTION 

Cosmic structure forms hierarchically, with smaller structures 
merging to form bigger ones. On the largest scales, clusters of 
galaxies merge releasing energies of the order of 10 64 erg on time 
scales of 1-2 Gyr (e.g. |Hoeft et aL]|2008| . During such merger 
events, large-scale shock waves with moderate Mach numbers are 
created in the intracluster medium (ICM). Shocks are collisionless 
features whose interactions in the hot plasma are mediated by elec¬ 
tromagnetic fields. They act on the ICM accelerating a fraction of 
the thermal distribution of particles transforming them into non- 
thermal populations of cosmic rays (CRs) through DSA (e.g. [En-| 
|sslin et al.|1998| ). 

Uncertainties remain on the mechanism required to pre¬ 
accelerate the electrons from thermal to supra-thermal energies 
before the injection in the DSA process, especially at low Mach 
numbers and in high beta plasmas (e.g. Kang et al. 2014, Vazza 


|& Bruggen||2014| ). Recent developments in the field have come 
from Particle-in-Cell (PIC) simulations that demonstrated that self¬ 
excitation of plasma waves via various instabilities may be crucial 
in the injection process jCaprioli & Spitkovsky|2014l|Guo, Sironi| 
|& Na rayan 2014]). A second mechanism that is expected to ac¬ 
celerate particles is turbulence amplified by cluster mergers (see 
Brunetti & Jones 2014 and reference therein). Evidence for the 
production of CRs in merging galaxy clusters by these two mecha¬ 
nisms has been found primarily in the form of radio relics and halos 
respectively. 


Radio relics are polarized, low surface brightness (^ 10 -6 Jy 
arcsec -2 at 1.4 GHz), steep-spectrum (a < -1, with S v oc y a ) 
sources that are typically located near the periphery of the cluster 
that underwent a recent merger. Relics are elongated and can reach 
sizes of 1-2 Mpc. The most likely scenario is that the synchrotron 
emitting electrons are (re)accelerated by outward moving shocks to 
a Lorentz factor of y e ~ 10 3 - 10 5 (e.g. 


Blandford & Eichler 


1987). 
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Radio halos are Mpc-size sources generated by radio-emitting 
plasma permeating the central parts of disturbed clusters. One of 
the potential mechanisms to form radio halos is the acceleration 
of relativistic particles by turbulence, generated in the ICM by 
cluster-cluster mergers (Brunetti et al. 2001; Petrosian 2001). Al¬ 
ternatively, the energetic electrons may be secondary products of 
proton-proton collisions (e.g. |Dennison|1980| . 

1.1 Double radio relics 

Models predict that giant radio relics typically occur along the clus¬ 
ter’s merging axis and their orientation should be perpendicular to 
the axis itself. Simulations of binary mergers ( |van Weeren et al.| 
2012b ) show that only mergers whose axes are fairly aligned with 
the plane of the sky can develop a double system of radio relics. A 
double radio relic system has a major advantage: the approximate 
knowledge of the merger axis orientation makes it easier to correct 
for projection effects and obtain the relic’s real distance from the 
cluster centre. Furthermore, when the merger axis is aligned with 
the plane of the sky, X-ray and optical observations are easier to 
interpret and cross-match with radio data. Only 15 other systems 
with double radio relics have been discovered so far |de Gasperin| 
|et al.|2014| . 

In this paper we report the discovery of a new double radio 
relic + radio halo system in PSZ1 G108.18-11.53 (z ~ 0.335). The 
Planck Collaboration (2014a) recently discovered this cluster by 
detecting its Sunyaev-Zel’dovich (SZ) signal. PSZ1 G108.18-11.53 
hosts one of the most powerful double radio relic system known 
to date and its relics have a very well defined morphology (non- 
disrupted arc-like structure, no strong blending with other sources, 
high surface brightness). These characteristics allowed us to per¬ 
form a detailed spectral analysis which led to the first observa¬ 
tional evidence of a non-uniform Mach number distribution along 
the shock that generated one of the relics. 

In Sec. [2] we present new radio and optical data. In Sec. [3] we 
present the global properties of the hosting cluster and of the radio 
relic (Sec. \3.2\ and halo (Sec. \3.3\ system. In Sec. [4] we analyse 
the radio spectra of the sources. Discussion and conclusions are in 
Sec. [5] and Sec. [6] respectively. Throughout the paper we adopt a 
fiducial ACDM cosmology with H 0 = 70 km s -1 Mpc -1 , Cl m = 0.3 
and Q a = 0.7. At the redshift of the target (z ~ 0.335) l"= 4.8 kpc. 
All regressions are made with an ordinary least squares algorithm 
that takes into account errors only in the dependant variable. Unless 
otherwise specified errors are at lcr. 


2 OBSERVATIONS 

The cluster was observed with the Westerbork Radio Synthesis 
Telescope (WSRT) at 21 cm (see Sec. IH) and subsequently with 
the Giant Meterwave Radio Telescope (GMRT) at three separate 
frequencies to enable spectral analysis (see Sec. |2.2| ). A summary 
of the radio observations is listed in Table [T] A preliminary opti¬ 
cal follow-up was also made using the W.M.Keck telescope (see 
Sec. [23}. 

2.1 Radio: WSRT 

The cluster was observed with the WSRT on January 2nd, 2014, 
for about 12 h using the default 21 cm set-up. Due to telescope 
upgrades, only 8 antennas participated in the observation. A total 
bandwidth of 160 MHz was recorded, spread over eight spectral 



Freq. 

(MHz) 

Obs. 

date 

Length 

(h) 

Resolution 

rms 

(juJy beam' 

GMRT 

147 

11 Jun 14 

8 

22"xl3" 

1600 

GMRT 

323 

02 Jun 14 

8 

8"x6" 

86 

GMRT 

607 

06 Jun 14 

8 

4"x4" 

43 

WSRT 

1380 

02 Jan 14 

12 

17"xl3" 

37 


Table 1 . Radio observations 

windows of 20 MHz in bandwidth. Each spectral window was fur¬ 
ther subdivided into 64 frequency channels. All four linear polar¬ 
ization products were recorded. The calibrators 3C286 and 3C147 
were observed at the start and end of the main observing run, re¬ 
spectively. 

The data were calibrated with CAS^] using the flux scale 
given by Perley & Butler ( 2013). The first step consisted of the re¬ 
moval of time ranges affected by shadowing. We then performed an 
initial bandpass correction using 3C147. Radio frequency interfer¬ 
ence was removed using the AOFlagger ( [Offringa, van de Gronde| 
|& Roerdink|2012} . This initial bandpass correction prevents flag¬ 
ging data affected by the strong bandpass roll-off. After flagging, 
we again calibrated the bandpass and subsequently obtained com¬ 
plex gain solutions for the two calibrator sources. The channel- 
dependent polarization leakage terms were calibrated using the un¬ 
polarized source 3C137. The polarization angles were set using 
3C286, carrying out a single correction per frequency channel. We 
then transferred all calibration solutions to the target source. For the 
target field, three rounds of phase-only self-calibration were per¬ 
formed. The data was imaged taking the spectral index into account 
during the deconvolution (i.e., nterms=2 |Rau & Cornwell|2011} . 
We used Briggs weighting with robust = 0 to gain sensitivity to¬ 
wards the diffuse emission. The final image was corrected for the 
primary beam attenuation. The final map is shown in Fig. [l]and has 
a resolution of 17"xl3" (rms: 37 //Jy beam -1 ). CASA is currently 
not able to perform a T sys correction; this could cause a few percent 
offset in the absolute flux scale. We cross-checked the flux values 
of 11 point sources with fluxes > 10 mJy visible in the NVSS sur¬ 
vey ( [Condon et al.|1998) finding our WSRT image systematically 
brighter by ~ 10%. We rescaled the final maps by this amount. 

2.2 Radio: GMRT 

For each observing frequency, the visibility data was processed us¬ 
ing an AIPS-basec0pipeline incorporating SPAM (Source Peeling 
and Atmospheric Modelling |Intema et al.||2009} ionospheric cal¬ 
ibration. The AlPS-based, semi-automated pipeline processing is 
very similar for each of the GMRT observing frequencies, being 
different only in the detailed settings of automated flagging. The 
pipeline uses the ParselTongue interface ( [Kettenis et al.||2006} to 
access AIPS tasks, files and tables from Python. Flux and band¬ 
pass calibrations were derived from 3C 286 after three iterations 
of flagging and calibrating based on the Scaife & Heald ( 2012) 
model. Additionally, instrumental phase calibrations were deter¬ 
mined by filtering out ionospheric contributions (see |Intema et al.| 
|2009} , an important step for direction-dependent ionospheric cal¬ 
ibration later on. Calibrations were transferred and applied to the 
target field data, simple clipping of spurious visibility amplitudes 
was applied, and data was averaged in time and frequency and 


1 version 4.2, http: //casa. nrao. edu 

2 AIPS: http://www.aips.nrao.edu 
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Figure 1 . Images of PSZ1 G108.18-11.53 at multiple frequencies obtained with WSRT and GMRT. Red contours trace the low-resolution version; halo and 
relics are labelled. Top-left: GMRT image at 147 MHz, black contours at (1,2,4, 8,16,32) x 3<x with cr - 1.6 mjy beam -1 (beam: 27"xl7"). Top-right: GMRT 
image at 323 MHz, black contours at (1,2,4, 8,16,32) x 3cr with cr = 86 /Jy beam -1 (beam: 1 l"x8"); red contour at 600 /Jy beam -1 (3<x of the low-resolution 
map, beam: 26"x20"). Bottom-left: GMRT image at 607 MHz, black contours at (1,2,4, 8,16,32) x 3<x with cr = 43 /dy beam -1 (beam: 6"x6"); red contour 
at 700 /dy beam -1 (3cr of the low-resolution map, beam: 30"x24"). Bottom-right: WSRT image at 1380 MHz, black contours at (1,2,4, 8 , 16,32) x 3<x with 
cr = 37 /dy beam -1 (beam: 17"xl3"). 


converted to Stokes I to speed up processing. The effective band- 
widths after flagging and averaging are 14.8, 31.3 and 29.2 MHz, 
centred on 147, 323 and 607 MHz, respectively. Self-calibration 
of the target field was started with an initial phase calibration us¬ 
ing a multi-point source model derived from the NVSS, WENSS 
and VLSS radio source catalogues ( Con don et al.|1 998 ; Rengelink 
|et al.|1997||Cohen et al.|2007| ), followed by (facet-based) wide-held 
imaging and CLEAN deconvolution of the primary beam area and 
several bright outlier sources (out to 5 primary beam radii). The 
visibility weighting scheme used during imaging is a mix between 
uniform and robust weighting (AIPS ROBUST -1), which gener¬ 
ally provides a well-behaved point spread function (without broad 
wings) by down weighting the very dense central UV-coverage of 
the GMRT. Self-calibration was repeated three more times, includ¬ 
ing amplitude calibration in the final round, and outlier Bagging on 


residual (source-model-subtracted) visibilities in between imaging 
and calibration. Next, two rounds of SPAM calibration and imaging 
were performed to minimize direction-dependent, residual phase 
calibration errors due to the ionosphere. In each round, (i) bright 
sources in and around the primary beam area were peeled (e.g., 
|Noordam|2004| , (ii) the peeling phase solutions were fitted with 
a time-variable two-layer phase screen model, (iii) the model was 
used to generate ionospheric phase correction tables for the grid of 
viewing directions defined by the facet centres, and (iv) the target 
held was re-imaged and deconvolved while applying the appropri¬ 
ate correction table per facet. We refer to |Intema et al.| < [2009]) for 
more details. 

Given the large range of baseline lengths of the GMRT for 
the 607 and 323 MHz observations we produced two maps, one at 
high resolution (robust = 0) and one applying a Gaussian weight 
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to the data (“tapering”). The final maps are shown in Fig. [T] For 
the 147 MHz map we obtained a beam of 27"xl7" (rms: 1.6 
mJy beam -1 ). From the 323 MHz observation we extracted an 
high-resolution map (beam: ll"x8", rms: 86 //Jybeam -1 ) and a 
low-resolution map (26"x20", rms: 200 //Jybeam -1 ) tapering at 
FWHM = 5 km. Finally, from the 607 MHz observation we ob¬ 
tained a high-resolution map (beam: 6"x6", rms: 43 //Jybeam -1 ) 
and a low-resolution map (30"x24", rms: 230 //Jybeam -1 ) taper¬ 
ing at FWHM = 2 km. 


2.3 Optical: Keck 

We conducted a spectroscopic survey of the target cluster with 
the DEIMOS instrument on the Keck II 10 m telescope on 2014, 
June 22. We used a similar instrument setup, slitmask design pro¬ 
cedure, and reduction procedure to that detailed in |Dawson et al.| 
( |2014| ). Observations were taken using 1" wide slits with the 
1200 line mm -1 grating, tilted to a central wavelength of 7400 A, 
resulting in a pixel scale of 0.33 A pixel -1 , a resolution of ~ 1 A 
(50kms -1 ), and typical wavelength coverage of 6000 A to 8600 A. 
The actual wavelength coverage may be shifted by ~ +410A de¬ 
pending where the slit is located along the width of the slitmask. 
For most cluster members this enabled us to observe Rj8, [OlII] 
4960 & 5008, Mgl (b), Fel, Nal (D), [Ol], Ha, and the [Nil] dou¬ 
blet. The position angle (PA) of each slit was chosen to lie between 
±5° to 30° of the slitmask PA to achieve optimal sky subtractioij^] 
during reduction with the DEEP2 version of the spec2d package 
( [Newman et al.|2013) . Within this range the slit PA was chosen to 
minimize the effects of chromatic dispersion by the atmosphere by 
aligning the slit, as much as possible, with the axis connecting the 
horizon, object and zenith (see e.g. |Filippenko|1982| . We observed 
a total of two slit masks with approximately 120 slits per mask. For 
each mask we took three 900 s exposures. 

Our primary objective for the spectroscopic survey was to 
maximize the number of spectroscopic redshifts of cluster mem¬ 
bers. This was challenging since the only available optical imag¬ 
ing data was from the Digitized Sky Surve}Q(DSS), which has a 
limiting magnitude of ~ 22 and ~ 3.5" full-width-half-max point 
spread function. Since the archival DSS catalogues do not push 
faint enough to detect a significant number of the cluster galax¬ 
ies we ran Sextractor ( |Bertin & Amouts||1996| ) on the DSS blue, 
red, and infrared digitized images following a similar setup and 
procedure to that detailed in |Jee et aL] ( [2015| . We found that the 
poor seeing of the DSS imaging coupled with the low galactic lat¬ 
itude (b = -11 deg) made star galaxy separation based on object 
shape impractical and selected spectroscopic targets based purely 
on color. We visually inspected the DSS color composite images 
to determine that blue - red > -2 appeared to exclude the major¬ 
ity of blue stars while being inclusive of apparent cluster galaxies. 
The slitmask creation procedure and reduction of the raw data are 
detailed in [Dawson et al.|(20l4) . 

We obtained 291 secure redshifts, however only 60 are of 
galaxies with the remaining 231 of stars; indicative of the impor¬ 
tance of deep subarc sec seeing imaging for spectroscopic target se¬ 
lection in low galactic fields. Of the galaxy spectra 42 are likely 
cluster members, defined as within 3000 km s -1 (three times the 
velocity dispersion) of the mean cluster redshift Cluster = 0.335. Of 


3 http://astro.berkeley.edu/~cooper/deep/spec2d/slitmask. 
html 

* https://archive.stsci.edu/cgi-bin/dss_form 



RA 02000) 


Figure 2. Composite image using the infrared, red and blue band of the 
DSS. Superimposed contours from the 323 MHz image obtained with the 
GMRT (contour levels as in Fig.[^ radio emission at 323 MHz). 


PSZ1 G108.18-11.53 


Redshift 

0 33466 +0 00061 

U.JJ4OO_ 000061 

Coordinates 

RA: 23 h 22 m 29.7 s - Dec: +48°46'30" 

X-ray luminosity a 

7.5 x 10 44 erg s” 1 

M 5 oo a 

(7.7 ± 0.6) x 10 14 M 0 

Velocity dispersion 

910+“ kms- 1 


a 0.1 - 2.4 keV, ROSAT 

b From SZ measurement, Planck Collaboration {2015} 


Table 2. Cluster properties 


the remaining spectroscopic galaxies 13 are foreground and 5 are 
background. 


3 PSZ1 G108.18-ll.53 


PSZ1 G108.18-11.53 was initially mentioned in [Planck Collabora¬ 
tion! ( |2014a| > but no redshift was known. In Planck Collaboration 


( 20 14b]) the cluster redshift was measured to be z = 0.336. Inde¬ 
pendently, with Keck observat ions , we measured a spectroscopic 
redshift of 0.33466^ ^} (Sec. 3.11. Throughout this paper we will 
use our redshift value. 

The cluster has been initially discovered by measuring its SZ 
signal over the cosmic microwave background radiation. The SZ 
signal integrated over the clusters angular extent measures the total 
thermal energy of the cluster gas, which is tightly related to cluster 
mass (e.g. Bor gani et al.|2 006). For PSZ1 G108. 18-11.53 this has 
been measured to be M 500 = (7.7 + 0.6) x 10 14 M 0 ( [Planck Collabo^] 
ration 2015). At these masses the Planck SZ-catalogue is complete 
at > 80% and only 72 more massive clusters are present in the cat¬ 
alogue. 

PSZ1 G108.18-11.53 is also detected in ROSAT (see Fig.[3j, 
which provide an approximate estimation of the total cluster lumi¬ 
nosity: L x ~ 7.5 x 10 44 erg s -1 . Unfortunately, the angular resolu¬ 
tion of ROSAT prevents us from carrying out a further analysis of 
morphological features. Using the (L, T) relations provided by [Pratt] 
|et al. | < f2009^ ) we estimated a cluster temperature of T~ 6.5 keV. 
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Counts per pixel 

0.1 0.2 0.3 0.4 



23m00.00s 40.00s 20.00s 22m00.00s 23h21m40.00s 

RA (J2000) 

Figure 3. Broad band (0.1 - 2.4 keV) photon image from the ROSAT all 
sky survey (contour levels as in Fig. [I] radio emission at 323 MHz). Image 
has been smoothed with a two-dimensional Gaussian kernel of cr = 2 pixels 
(1 pixel = 45 arcsec). Total estimated X-ray luminosity is Lx ~ 7.5 x 10 44 
erg s -1 . 


of 0.92. Thus the redshift distribution is consistent with being a 
Normal, providing no evidence of significant substructure along the 
line of sight. 

To examine the distribution of the cluster galaxies in projec¬ 
tion we use a color-magnitude selection based on the DSS photom¬ 
etry. We do not use the DEIMOS spectroscopic galaxies as there is 
a strong selection bias due to the limited number of spectroscopic 
masks observed and their limited field of view (5' x 16'). Instead 
we use the spectroscopic galaxies to inform a reasonable color- 
magnitude selection. Considering just the spectroscopic objects this 
selection results in 34 cluster galaxies, 5 foreground galaxies, 3 
background galaxies, and 67 stars. Thus there would be a large 
(Poisson) background component in any resulting galaxy density 
map, reflective of the poor imaging data. We show the photomet¬ 
ric galaxy number density map in Fig. [5] This map shows a gen¬ 
eral elongated over-density along the axis connecting the two radio 
relics, with the peak of the over-density located near the radio halo 
and the peak of the X-ray emission. This distribution is compatible 
with the idea of two sub-clusters which merged along a merger axis 
perpendicular to the line of sight and oriented along a line connect¬ 
ing the two radio relics. It is also interesting that the galaxy lumi¬ 
nosity density map shows a smaller over-density in the west, which 
may be related to the east-west elongated X-ray mission (Fig. [3j, 
although given the large background in the galaxy density map this 
is speculative. 



Figure 4. Distribution of spectroscopic measurements. Average value is at 
0 33466 +0 00061 

U. J J4UU_ 0 00061 . 

3.1 Galaxy distribution 

To investigate the distribution of the cluster galaxies we consider all 
42 galaxy spectra within the redshift range 0.3217 < z < 0.3483, 
which is Cluster ± 3 x cr, where z c iuster = 0.335 and cr is the ap¬ 
proximate velocity dispersion (lOOOkms -1 ). We show the redshift 
distribution of these galaxies in Fig. [4] We estimate the cluster red¬ 
shift and velocity dispersion using the biweight-statistic and bias- 
corrected 68% confidence limit ( [Beers, Flynn & Gebhardt||1990| 
applied to 100,000 bootstrap samples of the cluster spectroscopic 
redshifts. We find a cluster redshift of 0.33466^ ^} and a velocity 
dispersion of 910^° km s -1 . We caution the reader that care should 
be taken when interpreting these results given the small number 
of spectroscopic redshifts and that we have treated a likely major 
merger of two or more subclusters as a single system. We perform a 
Kolmogorov-Smirnov test to compare the redshift distribution rel¬ 
ative to the normal distribution defined by the calculated bi-weight 
location (redshift) and scale (velocity dispersion) and find a p-value 


3.2 Radio relics 

Two large radio relics are visible on the north-east and south-west 
side of the cluster. Their distance from the peak of the ROSAT X- 
ray emission is of ~ 1750 kpc and ~ 1280 kpc respectively. Their 
relative distance is of ~ 3100 kpc. The luminosity of the northern 
relic at 1400 MHz (^-corrected) is (26.6 + 0.1) x 10 24 W/Hz, while 
the southern relic luminosity is (18.0+0.1)xl0 24 W/H^] The north¬ 
ern relic is 33% more luminous than the southern although their ex¬ 
tension is similar (1500 and 1300 kpc, respectively). The northern 
relic brightness is not uniform across its extension peaking towards 
the south-east side. The external edge is slightly curved (a common 
feature for radio relics) while the internal edge is straighten Both 
internal and external edges show quite steep brightness gradients. 
The relic is thicker towards the south (390 kpc) and it becomes 
narrower towards north. Conversely, the southern relic brightness 
is uniform along its extension, but the source looks less compact 
and with less defined edges. The external edge has a quite steep 
brightness gradient, although half way along the relic extension, a 
filament of 240 kpc is visible in front of the relic. The relic thick¬ 
ness is uniform and of about 490 kpc. The surface brightness of the 
inner edge of the relic mildly decreases moving towards the cluster 
centre and its full extension is likely not entirely visible. 

A 10 - 30% polarization fraction is detected in the brightest 
regions both in the northern and southern relics (see Fig. [6}. Us¬ 
ing the map provided by Oppermann et al. ( 2012) we estimated a 
Faraday depth in the direction of PSZ1 G108.18-11.53 of -6.5 rad 
m“ 2 . The polarization angle has therefore been corrected for galac¬ 
tic Faraday rotation by adding an angle of 17.5°. The northern relic 
magnetic field appears aligned to the relic extension and curves 


5 Flux errors for extended sources are computed as S cn = cr x VW>eam> 
where cr is the local image rms and Abeam is the number of beams covering 
the source extension. 
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Figure 5. Cluster galaxy number density distribution (left) and unnormalized luminosity density distribution (right) in the field of PSZ1 G108.18-11.53 
(contour levels as in Fig. [I] radio emission at 323 MHz). Maps are based on DSS galaxies with a cluster membership selection trained on the spectroscopically 
confirmed cluster members. The map has been smoothed by a Gaussian with lcr radius = 50". White points are the cluster galaxies with spectroscopic data. 


Name 

Flux 147 MHi 
(mJy) 

Flux 323 MHz a 
(mJy) 

1 Flux 607 MHz a 
(mJy) 

Flux 1380 MHz a 
(mJy) 

Spectral 

index 

Luminosity 11 
(10 24 W/Hz) 

LLS 

(kpc) 

Relic N 

1032 ± 14 

422 + 2 

177 + 2 

67.7 + 0.3 

-1.25 + 0.02 

26.6 + 0.1 

1500 

Relic S 

753 ± 15 

323 ±3 

128 + 2 

45.4 + 0.3 

-1.28 + 0.02 

18.0 + 0.1 

1300 

Halo 

124+11 

56 + 2 

20+1 

6.8+ 0.2 

-1.40 + 0.07 

2.8+ 0.1 

850 


a fluxes are all extracted from a region equal to the 3 cr level of the 323 MHz map. 
b at 1400 MHz, k-corrected and spectral-index rescaled. 

Table 3. Radio properties 


close to the southern end. The southern relic polarization is less or¬ 
dered, but it shows again a magnetic field generally aligned with 
the relic extension. This magnetic field orientation might arise as a 
consequence of field compression by the shock wave which formed 
the radio relic (see e.g. [Iapichino & Brliggen|2012| >, it could come 
from magnetic field amplification from CRs turbulence fBruggen| 
|2013]>, o r it may be a prerequisite for efficient particle acceleration 
( |Guo, Sironi & Narayan|2014| ). 


3.3 Radio halo 

PSZ1 G108.18-11.53 has a centrally located radio halo. The source 
is quite uniform in brightness, and slightly elongated along the 
cluster merging axis. The halo maximum extension is ~ 850 kpc, 
with an average diameter of ~ 675 kpc|^] and a luminosity of 
(2.8 + 0.1) x 10 24 W/Hz at 1400 MHz (^-corrected). The halo has 
a rather steep spectral index of a = -1.40 + 0.07, close to the defi¬ 
nition of ultra-steep spectrum radio halo ( |Cassano et al.|2013) . The 
halo lies slightly below the cluster mass-power correlation found 
by |Cassano et al.| ( f2013j l as already noticed for other halos with a 
very steep spectrum. Given its luminosity, the radio halo is slightly 
smaller than expected from scaling-relations. [Cassano et aL]j2007| 
measured halo size and luminosity for 15 radio halos finding a clear 


Fractional polarization 

0.00 0.03 0.06 0.09 0.12 0.15 0.18 0.21 0.24 0.27 0.30 



Figure 6. Fractional polarization map with E-vectors displayed. The po¬ 
larization angle has been corrected for galactic Faraday rotation. Contour 
levels as in Fig.[^ radio emission at 1380 MHz. 


6 Measured following |Cassano et al.| {2007} as V^min x d max , d^ n an d 
dmax being the minimum and maximum diameter measured on the 3<x radio 
isophotes of the 323 MHz map. 


trend for more luminous radio halo to be more extended. For the 
luminosity of the halo in PSZ1 G108.18-11.53, the expected radio 
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Figure 9. Global spectral index analysis for the radio galaxies listed in Ta- 
ble[^(RG4 and RG6 are excluded). Solid lines are linear regressions, dotted 
lines connect data points. Fluxes have been arbitrarily rescaled to show all 
the spectra in a single plot retaining their shape in the log-log space. Error 
bars are at 5% of the data point value. 

halo average diameter is ~ 1000 kpc. Polarization is not detected 
down to the level of 10%. 

3.4 Radio environment 

The radio environment in the region of PSZ1 G108.18-11.53 is 
particularly rich with interesting sources. In Fig. [ 7 ] we show the 
full field of view of the GMRT 323 MHz observation. Apart from 
the two radio relics and the radio halo a number of extended radio 
galaxy (RG) are also present. The zoomed high and low-resolution 
images of these radio galaxies are in Fig. [8] while their properties 
are summarized in Table[4] Global spectral indices calculated using 
all four frequencies (147, 323, 607,1380 MHz) are shown in Fig. [9] 

RG1 Searching the NED archive we found a possible optical coun¬ 
terpart for RG1: 2MASX J23230751+4905386. Unfortunately the 
redshift of this galaxy is unknown but it is likely too far apart to 
be a cluster member and it is probably a foreground source. The 
global spectral index of this RG is a - -0.58 + 0.04 and its emis¬ 
sion is likely dominated by the central core. The extended lobes 
can be classified as an FR I with an interesting abrupt bending of 
one of the lobes similar to those of Z-shaped FR II radio galaxies 
( |Gopal-Krishna, Biermann & Wiita|2003] >. 

RG2 No optical counterpart has been found for this RG but given 
its extension and location it is likely a foreground object. This 
source is a rare example of a hybrid FR I/II radio galaxies with one 
FR I lobe and one FR II lobe (hybrid morphology radio sources, 
“HyMoRS” [Gopal-Krishna & Wiita|2000) . In depth analysis of this 
source will be presented in a forthcoming paper. 

RG3 No optical counterpart has been found for RG3 which is an 
FR II radio galaxy with a global spectral index of a - 0.77 + 0.09. 


The spectral index (see Fig. [9} shows signs of ageing as it becomes 
steeper at higher frequencies. 

RG4 Although no optical counterpart was found for RG4, its loca¬ 
tion and clear interaction with the ICM makes it a probable member 
of the cluster. The source is an head-tail radio galaxy whose lobes 
are bent by their interaction with the surrounding medium. 

RG5 RG5 has an identified optical counterpart located at 
RA:350.61°, DEC:48.69°(z = 0.333). The source is classified as 
a cluster member following the criteria of Sec. |2.3| The source is 
a wide angle tail radio galaxy, also common in the proximity of 
galaxy clusters. 

RG6 No optical counterpart has been identified for RG6. The 
source might be another example of HyMoRS although the sec¬ 
ond lobe does not show a clear hotspot presence. In depth analysis 
of this source will be presented in a forthcoming paper. 

RG7 RG7 has an identified optical counterpart in 2MASX 
J23212159+4827436. No redshift is available. The source appears 
to be a Z-shaped radio galaxy with prominent wings. Interestingly 
on the larger scales emission is also detected on the other side of 
the lobes compared to the wings direction. 


4 SPECTRAL ANALYSIS 

We extracted global spectral index values by fitting a first order 
polynomial to the integrated flux densities of the radio halo, the 
northern relic, and the southern relic (see Fig. ED- The radio halo 
spectral index value is a = - 1.40 + 0.07 but data points show a pos¬ 
sible curvature of its spectrum as expected by re-acceleration mod¬ 
els. The spectral index of the two radio relics is a = -1.25 + 0.02 
for the northern relic and a = -1.28 + 0.02 for the southern relic. 
Both spectral energy distributions appear flat, as seen for other ra¬ 
dio relics as long as the sampling does not span a wide enough 
range of frequencies (e.g. |Stroe et ak||2013) . However, as shown 
by Trasatti et al. (2015), integrated spectral index values of radio 
relics may be biased in a number a ways and their interpretation is 
not straightforward. 

The two-point spectral index map (Fig.|T0]) shows a clear gra¬ 
dient (flat —> steep) going towards the cluster centre along both relic 
widths. These gradients have been interpreted as a consequence of 
synchrotron ageing of electrons accelerated by a shock wave on the 
outer edge of the radio relic (e.g. |van Weeren et al.|2010| ), but such 
interpretation can mask more complex scenarios such as shocks 
with varying Mach number (Skill man et al.||2 013 ). While for the 
southern cluster, the steepening happens quite uniformly along the 
relic extensions, the northern relic’s outer edge has a different spec¬ 
tral index between the N and the S side, which may be related to its 
varying thickness and brightness along its extension. The spectral 
index of the radio halo is patchy, consistent with previous findings 
(Orru et al.|2007) . 

4.1 Plasma ageing 

Given the relatively high brightness of these radio relics we were 
able to perform a spectral ageing study of the two sources. Com¬ 
pared to radio galaxies, where this type of study has been made sev¬ 
eral times, only a limited number of radio relics have the required 
brightness. Moreover, in some cases the complex morphology of 
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Figure 7. Entire field of view at 323 MHz as seen by the GMRT, highlighted the regions zoomed in Fig. [8] Image noise at the field centre is 100 /rJy beam -1 , 
resolution is 10"x8". 


Name 

RA 

DEC 

Flux 147 MHz 

Flux 323 MHz 

Flux 607 MHz 

Flux 1380 MHz 

Spectral 


hh:mm:ss 

dd:mm:ss 

(mJy) 

(mJy) 

(mJy) 

(mJy) 

index 

RG1 

23:23:08 

49:05:39 

247 

174 

101 

67 

-0.61 ± 0.05 

RG2 

23:21:51 

48:50:23 

122 

77 

44 

13 

-0.99 + 0.17 

RG3 

23:21:24 

48:51:16 

106 

67 

39 

18 

-0.79 ± 0.06 

RG4 

23:22:27 

48:43:22 

- 

12 

6.5 

- 


RG5 

23:22:28 

48:41:25 

97 

44 

26 

13 

-0.91 ± 0.03 

RG6 (core/jet) 

23:21:46 

48:29:48 

< 16 a 

13 

9.8 

7.0 


RG6 (lobe N) 

23:21:48 

48:31:11 

- 

23 

17 

6.5 


RG6 (lobe S) 

23:21:36 

48:28:32 

- 

12 

> 7.4 b 

1.5 


RG7 

23:21:22 

48:27:43 

400 

232 

146 

76 

-0.74 ± 0.02 


a Source is blended with its northern lobe. 
b Source is only partially detected. 

Table 4. Radio galaxies 
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Figure 8. Zoomed images on some resolved sources in the 323 MHz field as observed with the GMRT. Given the different location across the map the rms 
values vary. Black contours trace the 607 MHz emission at (1,2,4, 8,16,32) x 3<x (from top-left <x = 62,47,47,62, 87 /Jy beam -1 , resolution 5"x5"), red 
contours show the emission at 323 MHz at (1,2,4) x 3<x (from top-left cr = 130,85,100,100,100 //Jy beam -1 , resolution 10"x8"). 
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Figure 10. Two-frequency spectral index map obtained between the two 
highest signal to noise ratio maps available: 323 and 1380 MHz. Contour 
levels as in Fig.[^ radio emission at 323 MHz. 



Figure 11. Global spectral index analysis for the two radio relic and the ra¬ 
dio halo. Solid lines are linear regressions. Error bars are at one <x calculated 
as explained in the text. 


the source made the spectral analysis of difficult interpretation (see 
e.g. |Owen et aT||2014[ for Abell 2256). Only recently extensive 
spectral age studies have been performed on the “Sausage” radio 


relic (Stroe et al. 

2013 2014aj) and on the “Toothbrush” radio relic 

(Ivan Weeren et a 

. 2012c). 


The flux calibration errors on all maps were set at 5%. To en¬ 
sure recovery of flux on the same spatial scales for the GMRT (147, 
323, 607 MHz) and WSRT (1380 MHz) observations we harmo¬ 
nized the uv -coverage and convolved the uniform-weighted maps 
to the same resolution of 22". We have used the BRATS package 


( [Harwood et al.|2013] > to fit the JP model ( |Jaffe & Perola|1973] > on 
a pixel-by-pixel basis. The JP model has been shown to be a good 
assumption for radio relics, or at least not worse than others ( |Stroe| 
|et al.|2013) . The radio images of the two relics were loaded sep¬ 
arately into BRATS, imposing a 3cr cut-off on the radio emission 
based on the off-source rms noise. The northern relic was covered 
by 196 pixels or an equivalent of almost 20 independent beams. 
The southern relic emission is fainter and part of the relic cannot 
match the cut-off criteria in shallower maps. The relic was covered 
with 115 pixels (~ 11 independent beams). 


The JP spectral ageing model has four free parameters: flux 
normalization, magnetic field B , injection index and time since 

acceleration f age . We assumed a uniform magnetic field of ~ 2 jdG 


(a reasonable assumptionjRyu et al.|20081 

Bonafede et al.|2010; van 

Weeren et al. 2010, Iapichino & Bruggen 

2012) which is able to re- 


produce the mass-luminosity scaling relation in double radio relics 


systems ( |de Gasperin et al.|2014| ). The shock aligns and likely en¬ 
hances the magnetic field in the vicinity of the shock front, but the 
magnetic field condition downstream of the shock front is largely 
uncertain and the assumption of uniformity is driven by our lack of 
knowledge. However, | van Weeren et al.| ( |2012c| ) showed that for the 
“Toothbrush” cluster magnetic field variations do not play a major 
role in creating spectral trends. Apart from synchrotron radiation, 
inverse Compton scattering of CMB photons contributes to elec¬ 
trons ageing. BRATS takes the inverse Compton effect into account 
as an additional magnetic field term B c mb added in quadrature to 
the original magnetic field B. The CMB energy density translated 
into B field terms is ( |Longair|2011| ): 

5 C mb = 3.18 (1 + z) 2 nG (1) 

where the redshift is z = 0.335 and the equivalent magnetic field 
for the inverse Compton interactions is B c mb = 5.7 jdG. 

The Mach number A4 of a shock is defined as 

Al — I’shock/ Csi (2) 

where c s is the sound speed in the upstream medium. For simple 
shocks, the a m j values can be used to infer the shock Mach number 
M following ( [Blandford & Eichler|1987| ): 

I 2ct in \ + 3 

M= ( 3 ) 


4.2 Injection index 


For each relic, we estimated an injection spectral index a m j by fit¬ 
ting the JP model with a range of a m j and leaving normalization 
and 4 ge free to vary. We made an initial injection index search be¬ 
tween 0.5 and 1.0 in steps of 0.1, followed by a fine search in steps 
of 0.01. The goodness of fit (sum of x 2 ) dependence on the as¬ 
sumed injection index has a single minimum (see Fig. [12] ), which 
represents the a m j that best describes the data. We note that it was 
common practice to estimate the injection index for a relic by look¬ 
ing at the flattest spectral index values close to the shock front or 
the integrated spectral index value. These approaches are subject 
to a number of limitations described in |Stroe et al.| < |2014a] ). Con¬ 
versely, estimating a m j as described here uses all the available data, 
minimizing the chances that a noise minimum/maximum in one of 

: the x 1 curve is 


12 


the maps can bias our result. As shown in Fig. 
well defined, hence the uncertainty in the injection index can be 
calculated following |Avni| (1976] ) who shows that the lcr error is 
given by a unit change in the x 1 value from its minimum. As we 
are considering the source on a pixel-by-pixel basis, the sum of x 2 
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Figure 12. Injection index aq n j for the entire northern relic (green circles), 
then N part of the northern relic (blue triangles), the S part of the northern 
relic (red squares) and the entire southern relic (cyan diamonds). Dotted 
lines are a 6th order polynomial fit to the data used to find the function 
minima. Bars in the lower part of the plot show the lcr error in the ct m] 
estimation. 


is over weighted by a factor of the beam area in pixel A beam . Thus, 
to find the errors, the sum of the has to be divided by A beam . An 
important caveat of this approach is the consequence of the limited 
resolution imposed by the beam in out maps. Mixing the emission 
from regions with different magnetic fields and/or different spectral 
ages, moves the spectra closer to power-law shapes, smearing the 
curvature ( |van Weeren et al.|2012c| ). 

As a next step, we fixed the injection value to the value we 
found: cr in j ;N = -1.06+gjJJ and to ar in j 5 s = -0.95+q^ for the north¬ 
ern and southern relic respectively. Finally, we fit a JP model with 
two degrees of freedom: normalization and f age . The resulting maps 
are in the first two columns of Fig. |T3] The southern relic dis¬ 
plays a quite uniform gradient from 0 to ~ 30 Myr across the 
fitted region. For the northern relic the situation is more com¬ 
plex: the region associated with a 4 ge = 0 (the shock front) ap¬ 
pears located only on the northern half of the relic while the aged 
electrons seem concentrated towards the south. This configuration 
would be hard to explain with current radio relic theories which 
point towards a shock present along all the relic extension, al¬ 
though we cannot exclude ad-hoc projection effects. However, as 
noticed in simulations, shocks do not have a uniform Mach num¬ 
ber along their extension ( [Skillman et al.|2013| ). Therefore, we at¬ 
tempted a separate fit for the a m j in the northern and southern half 
of the relic. The results are shown in Fig. [12] While the global fit 
found an cr in y N = 1.06^ ^, separate fits of the two regions returned 
a m j = 1.02+QQg and or in j = 1.17^^ for the northern and the south¬ 
ern half, respectively. The difference between the two injection in¬ 
dices is statistically significant. Fitting again for f age , this time fix¬ 
ing ofjnj to different values for the two relic regions, produced the 
map shown in the third column of Fig. |T3] The particle acceler¬ 
ation ( 4 ge = 0) happens now all along the external relic edge and 
the ageing happens along the relic shorter axis. It appears that some 
gradient is still present, but it would likely disappear if a continuum 
of cr in j was used. 

With the important assumption that the injection index - Mach 
number relation derived from DSA (Eq.[3j is correct, we estimated 
for the southern relic a Mach number of At = 2.33^^, while the 


northern relic Mach number goes from M ~ 2.20^ ^ of the north 
part down to M ~ 2.00^ ^ of the southern region, this can be 
the first observational evidence of varying Mach number along a 
cluster merger shock. 


5 DISCUSSION 


5.1 Injection index estimation 

An alternative approach to estimate the injection index is to use 
the average spectral index value of the shock region. However, 
this approach has some limitation as the definition of the region 
is arbitrary. In our case we used a region which covers the en¬ 
tire external edge of the relics, with a thickness of one beam. The 

spectral indices obtained in this way are: a in j ?N -1.1 + 0.1 and 

Oanj,s ~ -0.9 + 0.1. Meanwhile, using the integrated spectral index 
approach we would have underestimated the injection index, ob¬ 
taining a in y N = -0.75 + 0.01 and ar in j 5 s = -0.78 + 0.01. Many works 
in the literature use one of these two techniques to estimate ar in j 
(see e.g. Bonafede et al.|2012[ Giacintucci et al.|2008| Orru et al.| 


20071|Macario et al.|201 l[|Stroe et al.|2013||van Weeren et al.|2010| 


20l2c| for MACSJ1149.5+2223, MACSJ1752.0+4440, Abell 521, 
Abell 2744, Abell 754, CIZA J2242.8+5301 (“Sausage”), 1RXS 
J0603.3+4214 (“Toothbrush”)). In many cases, the inferred Mach 
number is higher than the Mach number measured from X-ray 
observations when they are available. For the Sausage relic, for 
instance, the temperature jump in X-ray observations indicates a 
Mach number of M ~ 3.2 ( |Akamatsu & Kawahara| 201 3]), lower 


than the radio Mach number of At ~ 4.6 derived by van Weeren 


letaklpOlOl ). On the other hand, for this cluster, by using the same 
method described in this paper, Stroe et al. ( 2014a) obtained a com¬ 
patible injection index a m j « 0.77 corresponding to M ~ 2.9. For 
the Toothbrush relic, the temperature and density discontinuities 
suggest M « 2 ( Ogrean et al. 2013), while from radio observa¬ 
tion the Mach number is Ad = 3.3 - 4.6. Recently |Hong, K ang 
& Ryu (2015) showed using cosmological hydrodynamic al simu¬ 
lations that this discrepancy might be due to an observational bias, 
as radio observations tend to pick up stronger shocks compared to 
X-ray data. Deep X-ray observation of PSZ1 G108.18-11.53 may 
unveil the presence of temperature discontinuity and will allow an 
independent estimation of the Mach number. 

In this paper we showed that with high-quality data we can 
find different injection indices along a single radio relic. Therefore, 
assuming the flatter spectral index region of the relic as representa¬ 
tive of the true injection index, would bias the estimation towards 
the region with the strongest Mach number (flatter spectra). Fur¬ 
thermore, a noise minimum/maximum in one of the maps might 
also bias the estimation towards flatter spectra ( [Stroe et al.|2014a| ). 
In an ideal case, the integrated spectra would be the best indicator 
as it would be marginally influenced by projection effects and uses 
all available data. However, the simple relation reintegrated = a inj +0.5 
does not hold in the case of evolving shocks ( jKang|20l5| ). Further¬ 
more, integrated spectral energy distributions of radio relics are not 
power laws when a large frequency range is considered ( |Trasatti| 
|et al. 12015] |Stroe et al.|2014b| >. This implies that the simple idea 
of a superimposition of spectra with an exponential cut-off (the so 
called continuous injection model) is not valid. Explanations for 
this include a non-constant magnetic field, a non-uniform Mach 
number, and/or that the electron acceleration process in low-Mach 
number shocks is not fully understood. 
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Figure 13. Particle age maps based on spectral ageing analysis, with uncertainty maps in the second column. Top and middle rows are for the northern and 
southern relic respectively; in the bottom row the injection index was fixed at two different values for the region north and south of the dashed lines (see 
Sec.fO). 


























The radio relic system in PSZ1 G108.18-11.53 13 


Surface brightness [mjy beam -1 ] 

3 6 9 12 15 18 21 24 27 30 



RA (J2000) 


Surface brightness [mjy beam x ] 

2.5 5.0 7.5 10.0 12.515.017.520.(22.25.0 



Figure 14. Comparison between the image at 323 MHz of MACSJ1752.0+4440 (left; |Bonafede et al.|2012) and of PSZ1 G108.18-11.53 (right). Contours 
in both cases are at (1,2,4,8,16,32) x 3cr with cr = 180 //Jy beam -1 and <x = 86 //Jy beam - Respectively. Red contour is the 3cr level for the low-resolution 
version of the maps. 


5.2 The sibling: MACS J1752.0+4440 

PSZ1 G108.18-11.53 appears from many points of view very sim¬ 
ilar to MACS J1752.0+4440 ( [Bonafede et al.|[20T2| |van Weeren| 
|et al.|2012a| ). Both clusters are at a similar redshift (z = 0.335 and 
z = 0.366), they are well-defined double radio relic + halo sys¬ 
tems, and their radio relics are the two most luminous ever discov¬ 
ered (see Fig. HD with a very similar luminosity. The two systems 
are also hosted in clusters with comparable SZ-derived masses: 
M s oo = (7.7 ± 0.6) x 10 14 M 0 and M 500 = (6.7 ± 0.4) x 10 14 M 0 , 
respectively. The two radio relic systems are also impressively sim¬ 
ilar from the morphological point of view. Both have a relic brighter 
than the other and they are not perfectly symmetric with respect to 
the cluster centre. Furthermore, both systems show similar degrees 
of polarization (~ 10%) with peaks of 30-40% on the two northern 
relics. The orientation of the magnetic field also looks compatible. 

A major difference is the relative distance between the relics. 
In the case of PSZ1 G108.18-11.53 this is ~ 3100 kpc, while for 
MACSJ1752.0+4440 it is around 2/3 of this, which may imply a 
more recent merger. Assuming an average shock Mach number of 
2 for both clusters as derived from radio observations, we can infer 
the upstream temperature from the Rankine-Hugoniot relations for 
a T = 5/3 gas to be 



Here T u is the upstream temperature and T d is the downstream tem¬ 
perature which we equate to the mean cluster temperature. The 
upstream temperature can then be converted into a sound speed 
c s = 1480 yjT /(10 8 K) kms -1 , yielding c s ~ 870 km s -1 for PSZ1 
G108.18-11.53 and c, ~ 880 kms -1 for MACS J1752.0+4440. Fi¬ 
nally, combining this information with the relic distance we can cal¬ 
culate the collision tim^Jto be ~ 870 Myr ago for PSZ1 G108.18- 
11.53 and ~ 530 Myr for MACS J1752.0+4440. Likely the mean 

7 With “collision time” we mean the time since the two cluster cores have 
reached their minimum separation. This time is « 1 Gyr shorter than the 


cluster temperature is lower than the post-shock temperature (e.g. 
|Markevitch|2006) , so these ages should be considered upper limits. 

The two relics in PSZ1 G108.18-11.53 are similar in di¬ 
mensions, while the relics in MACS J1752.0+4440 are asymmet¬ 
ric in power with the flux ratio between the two relics being 
larger in MACS J1752.0+4440. The most striking difference be¬ 
tween the two system is the radio halo which is centrally located 
and disconnected from the relics in PSZ1 G108.18-11.53, while 
in MACS J1752.0+4440 it extends touching both radio relics ([van] 
Weeren et al. 2012a). Possible explanations include: 

(i) The radio halo in MACS J1752.0+4440 is at its maximum 
extension while the halo in PSZ1 G108.18-11.53 is already in its 
fading phase, in line with its relatively steep spectral index ( |Don-| 
Inert et al.|2013) . 

(ii) A different merger dynamics caused different amounts of 
turbulent energy to be injected into the ICM. 

(iii) The radio halo follows the magnetic field which traces 
the X-ray distribution. The X-ray brightness distribution appears 
more compact in PSZ1 G108.18-11.53 while it is elongated in 
MACS J1752.0+4440 jvan Weeren et al.|2012aj ), however this will 
need to be confirmed by Chandra observations. 

(iv) The spatial distribution of seed electrons could differ in the 
two clusters. We do not know yet how CRs accumulate in clusters 
over time, but they live for ~ 1 Gyr and at very low levels (0.1 %) are 
expected to be injected from supemova/AGN feedback ( [Brunetti &| 
|Jones|2014) . Hence, this explanation only shift the problem towards 
the seed particles. 

If indeed radio relics behave self-similarly, i.e. similar merger 
configurations produce similar radio relics, the origin of radio relics 
should be generic. If radio relics require shock re-acceleration (e.g. 
|Bonafede et al.|2014) , the pre-existing radio dark CR populations 
also needs to be similar across clusters. This would make it unlikely 


time since the centre of mass of the subcluster crosses the virial radius of 
the main cluster |Tormen, Moscardini & Yoshida|2004| . 
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Figure 15. Cluster mass (M 500 ) versus the hosted double radio relic sys¬ 
tem luminosity. The correlation in this figure is discussed in |de Gasperin| 
|et al. 1|20 14) , the red circle is PSZ1 G108.18-11.53 while the blue circle is 
MACSJ 1752.0+4440. 


that the fortuitous co-location with a source of seed CR -as an aged 
radio galaxy lobe- is a requirement for the presence of a radio relic. 


6 CONCLUSIONS 

We reported the discovery of a new double radio relic + radio halo 
system in PSZ1 G108.18-11.53. The system is one of the (now) 16 
systems with two radio relics. Among these systems only six others 
also show a clear radio halo. Our main results are summarized here: 

• PSZ1 G108.18-11.53 hosts the second most powerful radio 
relic system known, with a total power of Li. 4GHz = (44.6 + 0.3) x 
10 24 W/Hz. The cluster also hosts a radio halo of power Li. 4GHz = 
(2.8 + 0.1) x 10 24 W/Hz. Both, radio relic and halo emission agree 
with luminosity-mass scaling relations found by ( jde Gasperin et al.| 
|2014| > and ( [Cassano et al.|2013) , respectively. 

• We measured up to ~ 10 - 30% polarized emission from both 
the northern and southern relic. The magnetic field in the northern 
relic appears well aligned with the relic extension. 

• We measured the injection index of the radio spectrum for 
both relics. We found that the northern relic local spectra are bet¬ 
ter described by two different injection indices, a m j = 1 . 02 +gg 4 and 
crinj = 1 . 17+0 07 f° r the northern and the southern half, respectively. 

• The Mach number derived from the injection index accord¬ 
ing to simple DSA theory (see Eq. [3j has been measured to be 
At = 2.33+026 f° r the southern relic, while the northern relic 
Mach number goes from At = 2.20 +q^ in the north down to 
M = 2 . 00 !»;“ in the southern region. This indicates a gradient 
in the Mach number along the shock which produced the northern 
relic, as predicted by simulations (see e.g. ISkillman et al.|2013) . 
The Mach number across the relic varies by ~ 10%. 

• We note the interesting similarity between 
MACS J1752.0+4440 and PSZ1 G108.18-11.53, the first and 
second most luminous radio relic systems known. They are 
at similar redshift, have similar host masses, and have similar 
radio luminosity and relic relative position and morphology. An 
important difference between the two systems is the relic distance, 
MACS J1752.0+4440 being 2/3 of PSZ1 G108.18-11.53, and the 
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Figure Al. Total intensity map obtained using the MSMF algorithm as de¬ 
scribed in Sec. [A] The reference frequency is 380 MHz and beam size is 
12"xll". Contours are at (1,2,4, 8,16,32)x3cr with cr = 100 /Jy beam -1 . 


radio halo which connects the two relics in MACS J1752.0+4440 
( [Bonafede et al.||20l2| |van Weeren et al.||2012a| ), while in PSZ1 
G108.18-11.53 is confined in the central part of the cluster. 


APPENDIX A: GMRT IMAGING USING MSMF 
ALGORITHM 

The novel multi-scale multi-frequency deconvolution algorithm in¬ 
cluded in latest CASA releases, models the wide-band sky bright¬ 
ness distribution as a combination of both spatial and spectral func¬ 
tions { [Rau & Cornwell|2011| ). The algorithms is therefore capable 
to perform radio image reconstruction on wide bandwidth obser¬ 
vations or on observations carried out at multiple frequencies. The 
output of the process is a monochromatic brightness distribution 
map at a chosen frequency and a map of spectral indices extracted 
from the Taylor series used to model the spectra of the flux com¬ 
ponents. We jointly deconvolved three GMRT datasets at 147, 323 
and 607 MHz, modelling the source spectra with a fist order poly¬ 
nomial. To our knowledge this is the first time this approach is used 
on GMRT data. 

We show the outcome in Fig. m where we produced an im¬ 
age at the intermediate frequency of 380 MHz. The image has a 
comparable resolution with the 323 MHz image (Fig. [TJ. Despite 
the slightly higher noise, the extended emission appears better re¬ 
constructed due to a better coverage of the uv- plane. We also ob¬ 
tained the spectral index map (Fig. |A2| ) which appears qualitatively 
similar to the 2-point spectral index map shown in Fig. [lO] The 
2-point spectral index map was computed using the two highest 
signal-to-noise ratio maps (323 and 1380 MHz). It is not surpris¬ 
ing that sampling higher frequencies than in the MSMF case, the 
spectral index values are on average ~ 0.1 steeper for the relic and 
~ 0.2 for the halo which has a more curved spectrum. 
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Figure A2. Spectral index map obtained using the MSMF algorithm as de¬ 
scribed in Sec.lA] Blanked pixels had an error higher than 0.5. Contour 
levels as in Fig.^radio emission at 323 MHz. 
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